Multiple Time Scales in Diffraction Measurements of Diffusive Surface Relaxation 
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We grew SrTi03 on SrTiC>3 (001) by pulsed laser deposition, using x-ray scattering to monitor 
the growth in real time. The time-resolved small angle scattering exhibits a well-defined length 
scale associated with the spacing between unit cell high surface features. This length scale imposes 
a discrete spectrum of Fourier components and rate constants upon the diffusion equation solution, 
evident in multiple exponential relaxation of the "anti-Bragg" diffracted intensity. An Arrhenius 
analysis of measured rate constants confirms that they originate from a single activation energy. 

PACS numbers: 61.10.-i,68.55.Ac,81.15.Fg 



Despite decades of study, a complete description of 
the fundamental mechanisms that control a broad class 
of epitaxial growth techniques persists as a challenging 
problem in non-equilibrium physics. For example, in the 
thin film-growth technique of Pulsed Laser Deposition 
(PLD), an extremely dense plume of laser-ablated ma- 
terial deposits a high concentration of ionic and neutral 
species onto an atomically smooth surface Numerous 
excellent diffraction-based studies @> 0, 0> U°f the ensu- 
ing kinetics have revealed multiple time scales in the post- 
pulse intensity transient, and have suggested that two or 
more distinct physical mechanisms govern the surface's 
relaxation toward equilibrium. In this Letter, we combine 
elementary diffusion theory and kinematic scattering the- 
ory with experiment to show that a single energy process 
can produce a non-exponential diffraction response dur- 
ing PLD. A better understanding of surface kinetics may 
explain why PLD can produce atomically smooth films, 
rivalling other growth techniques (e.g. molecular-beam 
epitaxy, chemical vapor deposition). The sensitivity of 
functional properties to smoothness necessitates this un- 
derstanding, since PLD is a candidate to grow device- 
quality films of complex-oxide materials 0, 0, @ • 

Film-growth models typically assume that particles 
diffuse on the substrate until they either evaporate, at- 
tach to an existing step edge, or run into other particles 
and nucleate islands [§, . Depending on the substrate 
temperature and specific energy barriers, films grow in 
one of several well-known modes: 3D, layer- by-layer, or 
step- flow 0, 0, 0| . These broad categories encompass 
myriad film morphologies, and the specific morphology 
determines the form of the intensity in diffraction stud- 
ies of film growth. 

To obtain information regarding morphological evo- 
lution, we measured the Small-Angle X-Ray Scattering 
(SAXS) during PLD (Fig.0}. A CCD detector recorded 
both the x-rays reflected specularly from the surface, and 
those with a small in-plane component, scattered in a 
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FIG. 1: X-ray intensity measured around (00|) during 
SrTi0 3 homoepitaxy at T = 970 K, Po 2 = 10" 5 Torr and 
/ = 0.1 Hz. An absorber attenuated the specular beam. In- 
set: Diffuse intensity during deposition of first half-monolayer. 
Each curve represents the integrated signal following a pulse. 



narrow angular range about the specular beam. Because 
the specular intensity dwarfed the diffuse by three or- 
ders of magnitude, we attenuated the specular with an 
aluminum absorber so that both signals could be simulta- 
neously measured without saturating the CCD detector. 

To increase the diffuse signal, we measured the x-rays 
scattered not around the often-used "anti-Bragg" (00^) 
position, but rather around the (00-i) position of the 
specular Crystal Truncation Rod (CTR). In general, the 
specular and diffuse signals intensify as one moves along 
the CTR toward (000), although the intensity profile de- 
pends strongly upon the surface morphology In 
layer-by-layer homoepitaxy, the specular intensity along 
the entire CTR oscillates once as each new layer is grown. 
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During layer-by-layer SrTi03 (STO) homoepitaxy, we 
observed that the diffuse and specular signals oscillate 
at the same frequency, out of phase (Fig. P). During 
growth of the first monolayer, the diffuse peaks are sep- 
arated by A(j|| « 0.06 A" 1 (Inset to Fig. [TJ. This indi- 
cates that surface features with a characteristic length 
scale of ~ 200 A develop as the growing layer nears half- 
completion, and vanish at layer completion. This length 
scale increases with increasing film thickness, eventually 
preventing resolution of the diffuse and specular signals. 

We interpret these data to indicate the nucleation 
of unit-cell high islands separated by ~ 50 unit cells 
(a-STO = 3.905 A). This morphology is consistent with 
20 nm round features visible in our atomic force micro- 
scope ( AFM) images and those in the literature 0] . The 
nearly constant value of q» (during the first monolayer) 
of the diffuse peaks indicates that while the islands grow 
in size as new material arrives, little subsequent nucle- 
ation occurs after the first pulse. As islands coalesce, the 
diffuse signal corresponds to holes in the surface that fill 
in as the specular intensity approaches a maximum. 

From the specular oscillations of Fig. ^ we observe 
that growth proceeds at roughly one monolayer (ML) 
per 100 s. The laser repetition rate (/ = 0.1Hz) and cov- 
erage of deposited material (ctq wO.l ) set this time 
scale. The shorter time scales associated with material 
incorporation and diffusion are regulated by the substrate 
temperature T, surface morphology, and diffusion energy 
barrier U (and possibly also by the kinetic energy of inci- 
dent particles 0|). While a rate-equation model 0,0] 
provides much qualitative insight regarding diffusion, the 
following continuum model allows direct extraction of U. 

In a continuum model of diffusion, the surface density 
p of diffusing species evolves according to the diffusion 
equation: || = DV 2 p . The diffusion coefficient, D, de- 
pends on T and U, via D oc e ~ u / kBT 0. The Fourier 
transform of the diffusion equation yields p q cx e~ Dq 
where q corresponds to the spatial frequencies present 
in p. In a boundary value problem, one expects a dis- 
crete spectrum of rate constants koc—Dq 2 . Short wave- 
length modes decay rapidly, and the spatial geometry 
determines the fc-spectrum. We emphasize that a single 
energy U gives rise to multiple diffusion rates k. 

The combined SAXS (Fig. [TJ and AFM data suggest 
a film morphology of round islands and holes. We there- 
fore solve the diffusion equation in a circular region of the 
(non-miscut) substrate of radius rt,, with a step located at 
r a < Tb- We write the solution as p(r, c/>, t) = X(r, <f))er kt , 
and specify X < and X> as the spatial solutions in the 
regions < r < r a and r a < r < r^, respectively. As- 
suming || = 0, X takes the form of Bessel, Jo(r), and 
Neumann, No(r), functions. 

We assume the step acts as aperfect sink, providing 
a boundary condition p(r a ) = 0,llfj. This method ne- 
glects the possible energy barrier for particles to cross 



downward over the step. One can incorporate this 
"Ehrlich-Schwoebel barrier" by speci fying the particle 
current j — —DVp at the boundary |l8L llflj. without 
changing the conclusions of this paper. 

The X < exclude the AT (r), which diverge as r^0: 

A< = A m J (a m r) = ^A m< . 

m tci 

Here a m r a is the mth root of Jo{r). 

We further assume that §j?|r& = 0, i.e. that the system 
boundary is far from any steps. Then 

A> = B m Mj3 m r) + C m N o [0 m r) = ^ x ™> , 



where the boundary conditions determine the (3 m and the 
ratio jf^. Finally, we specify p(t — 0,r) — <To, where do 
is the coverage of single pulse. This sets A m and B m via 
the orthogonality of the Jo(r) and No(r). 

Kinematic theory accurately describes the scattered in- 
tensity. For the simple case of momentum transfer nor- 
mal to the surface of a non-miscut crystal, one has, in 
homoepitaxy |2fj|. 
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Here q is the scattering vector, F (q) is the scattering 
amplitude of a single layer, 9 n (t) is the time-dependent 
coverage of the n th layer, and d is the layer spacing. 
The first term represents the scattering from an ide- 
ally terminated single crystal. The second term rep- 
resents the scattering from the deposited film. At the 
anti-Bragg position, x-rays scattered from adjacent lay- 
ers interfere destructively, providing maximum sensitiv- 
ity to single step height fluctuations. Researchers often 
exploit this property in reflection high-energy electron 
diffraction (RHEED) 0, S S E3 and x-ray growth stud- 
ies 

Combining the complete solution p(r,t) with Eq. QJ, 
we calculate the anti-Bragg x-ray intensity, J(t): 



I(t) <x(^-8 -a 



2]>> m e-^ 



(2) 



where Co is the coverage of single pulse. The pre-pulse 
surface coverage, 6*o, the coverages a m of species in the 
mass-losing layer, and the choice q m = a m or {3 m de- 
pend on whether the step at r a bounds an island or a 
hole. The complicated time dependence of Eq. J2J in- 
volves multiple rate constants k m = Dq^, rather than a 
single exponential. 

At the nucleation of each new layer, the morphology 
is characterized by islands. In this regime, 9o = (t^) 2 j 
c ot = %i J ' a X m< (r)rdr, and q m = a m . Only species 
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FIG. 2: Calculated anti-Bragg x-ray intensity following a 
laser pulse, using Eq. @ and parameters in agreement with 
experiment. The normalized intensity drops instantaneously 
upon pulse arrival, and then either falls (a) or rises (b), de- 
pending upon surface morphology. Tables in each case list 
fractions of total deposited mass, ^ , and diffusion rates of 
the first four modes. Integral forms of a m were evaluated 
numerically. 
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FIG. 3: (a) Fits of Eq. © to 900 K data with one (dashed 
curve) and three (solid curve) rate constants. The circle size 
roughly equals the error of each point, (b) Residuals of the fit 
with three rate constants. Inset: Arrhenius plot of 13 samples, 
with T = 900 -> 1060 K, and Po 2 = 10" 6 Torr. Ordinate is 
ki (diamonds) and (squares) of Eq. ©. 



transferring off of the island affect /, so that X> does 
not influence the time-dependence. 

In this case, / slightly decreases following the large in- 
stantaneous drop due to pulse arrival [Fig. Ufa)]. The 
growing layer moves toward half-completion and a min- 
imum in the anti-Bragg oscillation. The diffusion- 
moderated change is small, since most of the deposited 
material does not land on the island. The small fractions 
of total deposited material, ^fN transferring between lay- 
ers, move at relatively rapid rates k m [Inset to Fig.^a)]. 
Many researchers have observed intensity transients in 
qualitative agreement with this form 0, . 

As the islands grow and coalesce, a network of holes 
better describes the surface. Here, 9 = 1 — (^J-) 2 , 
f m = pr JJ? X m> (r)rdr, and q m = f3 m . The intensity 
rises substantially following the instantaneous drop, as 
the film moves closer to the layer completion condition 
and growth oscillation maximum [Fig.^b)]. The k m are 
in general slower than in the island case. Roughly 80% of 
deposited material is in the mode associated with X 0> , 
with a time constant equal to a few seconds [Inset 
to Fig. Ufa)]. Of the species transferring into the hole, 
~ 10% move at rates k m ;§> fco, adding faster components 
to the intensity change, noticeable at early times. 

To test the model, we measured the anti-Bragg in- 
tensity during STO homoepitaxy via PLD (Fig. [3J, un- 
der conditions which promote steady-state layer-by-layer 
growth [l3j | . Our PLD chamber is an integral component 
of a fully featured x-ray diffractometer that is perma- 
nently installed in the G3 hutch of the Cornell High En- 
ergy Synchrotron Source (CHESS) jig. By controlling 



the chamber oxygen pressure, Po 2 i substrate tempera- 
ture, T, and laser repetition rate, /, we can select differ- 
ent growth modes, while monitoring film growth in-situ 
with x-rays. We grew STO films on STO (001) substrates 
with P 02 = 10~ 6 Torr, T= 900 -> 1060 K, and / = 0.1Hz. 
We set the laser energy density at the single crystal STO 
target to 2 — 3 , and the target-substrate distance to 
6 cm, for an average er « 0.1 ^jj^. To prepare a Ti0 2 - 
terminated surface, we etched the substrates in buffered 
NH 4 F-HF [H|23. We annealed each sample in 1 mTorr 
O2 for 1-2 hours at 1060 K prior to growth to obtain a 
smooth surface |24j . 

To obtain the required time resolution, we recorded 
I(t) with a multi-channel scaler (MCS) set to 10 ms dwell 
time per channel. The MCS triggered the laser at the 
midpoint of each pass, collecting 5 s of intensity before 
and after each laser pulse. We then normalized each MCS 
pass by I(t) integrated over the pre-pulse region |l6| . 

We typically measured anti-Bragg intensities of 10 4 ^p. 
This signal strength yields 10% statistical noise in each 
MCS bin, comparable to the signal change ^j-. To reduce 
the noise, we averaged normalized MCS passes for a given 
sample over the entire deposition jlfj. While improving 
counting statistics, this method sacrifices information by 
combining intensity transients from a range of surface 
morphologies. For example, during the growth of the 
first half-monolayer, individual transients resemble the 
curve of Fig. Ota) (noise prevents quantitative analysis). 

However, sums of MCS passes at constant growth os- 
cillation phase reveal the transient to resemble the curve 
of Fig. 12ft)) at all phases 0] . This suggests that a hole 
morphology dominates throughout much of the growth. 
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A numerical simulation of a system containing an island 
and a hole reveals that the holes dictate the time struc- 
ture of the intensity transient. One can crudely predict 
this result by summing the two curves of Fig. |3 Nu- 
merous authors have observed rising intensity transients 
during the "falling" portion of the anti-Bragg oscillation, 
where one might expect an island morphology 

To fit the data, we truncated the series of Eq. (J2J) at 
m = 3, for which the model predicts the fe m >3 to be much 
greater than /Nyquist = 50 Hz. In the hole regime, we 
have a m <C 1. Retaining small quantities to first order, 

3 

/(*) = a - J2 bme- k ^ . (3) 

m— 1 

We fit Eq. to data from 13 samples grown un- 
der layer-by-layer conditions. For samples grown at 
T>970K, we required only two rate constants, suggest- 
ing that k 3 > /Nyquist at high T. A typical fit [Fig. Eta)] 
describes the data well, with xl = 1-1 and normally dis- 
tributed residuals [Fig.Efb)]. Adding terms m> 3 did not 
improve xt- For comparison, we also show a fit with the 
series truncated at m = 1, a form often assumed for the 
intensity evolution 0, Q . This simple exponential clearly 
fails to describe the data at early times (xl = 4.9). 

Averaged over all samples, the ratio |^ = 11.9 agrees 
with the hole model prediction of 9 to 13, which can vary 
depending on hole size and shape. The ratio |^ = 80 
exceeds the hole model prediction of 30. This may be 
due to a small island contribution, for which fc , 2 '"' a "' i ss 
100. Additionally, the fixed r a approximation used to 
derive Eq. J2J decreases the ratio n ^ >1 ■ In reality, the 
holes fill in, effecting longer saturation times for the x-ray 
intensity, and reducing the measured k\. 

The ratio || s=s 2 is smaller than the model prediction 
of 10. High frequency terms thus contribute more to 
the signal than predicted. This discrepancy could arise 
from the approximation of a two-level system. The actual 
broad growth interface would place more material near 
step-edges, increasing the rate of the intensity change. 

In the inset to Fig. EJ we present k\ and ki in an Ar- 
rhenius plot for the 13 samples mentioned above. Since 
k m oc D, the slopes of linear fits give the diffusion activa- 
tion energy, U. The slopes of separate fits to the two dis- 
tributions agree to within 1%, revealing that both rates 
arise from a single energy. A simultaneous fit of both 
distributions yields U = 0.6 ± 0.2 eV, in the range of 
published diffusion barriers for oxides |2fij . One source 
of error, due to sums over multiple MCS passes could be 
eliminated by higher x-ray intensity, allowing analysis of 
individual laser pulses. 

That a single process generates multiple time scales 
in the diffracted intensity is a direct consequence of the 
characteristic length scale on the surface, revealed in our 
SAXS data. The presence of this length scale forces a 
discrete q-spectrum upon the density of diffusing species. 



Higher-q components decay quickly as particles diffuse to 
steps, effecting a rapid early response in the time-resolved 
anti-Bragg intensity. This effect is more easily seen in 
PLD than in steady state growth techniques, due to the 
high density immediately following a laser pulse. 
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